Signatures of a gap in the conductivity of graphene 
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The electrical conductivity of suspended graphene has recently been measured for the first time, 
and found to behave as a ~ ^/\n\ as expected for Dirac quasiparticles at large carrier density. 
The charge inhomogeneity is strongly reduced in suspended samples, which revealed an unexpected 
insulating trend in o(T). Above a transitional density n* , the temperature dependence was found 
to revert to metallic. We show that these features of the DC conductivity are consistent with a 
simple model of gapped Dirac quasiparticles, with specific signatures of a gap in the vicinity of the 
charge-neutrality point. Our findings are reminiscent of the conductivity profile in semiconducting 
materials, exhibiting a thermal activation for T > T and a weakly T-dependent background for T < 
T, where T is given by the saturation density h associated with the residual charge inhomogeneity. 
We discuss possible origins of a bandgap in graphene as well as alternative scenarios. 
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PACS numbers: 73.63.Bd, 71.30.+h, 05.10.Ln 



The experimental study of graphene has recently pro- 
gressed towards an improved understanding of the intrin- 
sic properties of this new carbon nanomaterial. Notably, 
a clear signal for the fractional quantum Hall effect has 
been found [l|, |2| , and carrier mobilities far in excess of 
silicon-based devices have been reached [3J. While the 
high quality of graphene and the Dirac nature of the 
charge carriers is by now well established, the issue of 
gap formation has remained far less clear. Indeed, sev- 
eral features suggestive of a gap but difficult to interpret 
have been reported under various circumstances. The 
Lanzara group has performed ARPES studies [J, [5| of 
graphene layers on a substrate, whereby a gap-like fea- 
ture was detected and attributed to substrate effects. On 
the other hand, the conductivity of suspended graphene 
devices was characterized by the Kim group [6| and the 
Andrei group [7] . In addition to demonstrating a strongly 
T-dependent DC conductivity er, a transitional carrier 
density n* was found, where a(T) changes from metal- 
lic to insulating. A different perspective is provided by 
the Andrei group in Ref. Q via STM spectroscopy of de- 
coupled graphene flakes on graphite, where a ~ 10 meV 
gap centered on the Dirac point was found at zero mag- 
netic field, accompanied by a corresponding splitting of 
the lowest Landau level. 

On the theoretical side, a quantitative explanation 
of a(n,T) has remained elusive. For instance, the ap- 
pearance of a transitional density n* is incompatible 
with electron-phonon scattering [y, |9(, while electron- 
electron scattering leads to metallic behavior at the neu- 
tral point [lfj . A qualitative description of the empirical 
cr(n, T) is give n by Landauer transport theory for ballistic 
graphene [11| . The Hall probe lead geometry of Ref. [6[ 
minimizes the effects of a finite sample size, whereas the 
two-lead geometry of Ref. [7| is better described in the 
Landauer approach. Our primary objectives are to deter- 
mine whether the measured cr(T) of Ref. [fjj at the neu- 
tral point can be quantitatively described in terms of free 



gapped Dirac quasiparticles and a weakly T-dependent 
background, and whether such a description is consistent 
with <r(n, T) at finite carrier density. 

The rationale for the gapped quasiparticle scenario is 
as follows: Initially, the suspended samples of Ref. [fj 
showed an essentially T-independent DC conductivity, 
which upon current-annealing acquired a pronounced T- 
dependence of insulating type. While the increase in re- 
sistivity at low T remained modest (a factor of ~ 3 in the 
range 5 — 150 K), the data suggests the existence of two 
regimes, reminiscent of conventional semiconductors [12j : 
a low-T regime (5 — 35 K) where a is dominated by a 
weakly T-dependent background a, , and a thermally ac- 
tivated regime (35—150 K) where a increases rapidly with 
T. Notably, several mechanisms for gap generation have 
been proposed, such as explicit breaking of the sublat- 
tice symmetry by mechanical strain [13j . or spontaneous 
induction of a Mott insulating state via strong Coulomb 
interactions [14| . It is plausible that a b is due to residual 
charge density inhomogeneities [15] which, however, are 
reduced by an order of magnitude compared with sam- 
ples on a substrate [7]]. Regardless of the origin of a b , 
we find that it can be reliably subtracted. 

The Hamiltonian describing Dirac quasiparticles with 
a gap A and Fermi velocity v F ~ c/300 is given by 



H = a 1 v F k 1 + a 2 v F k 2 + <r 3 A/2, 



(1) 



where the cr, are Pauli matrices, and we account for the 
spin and valley degeneracy below. From this starting 
point, a is calculated as the diagonal part of the conduc- 
tivity tensor c , given by the Kubo formula 
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where j3 = \/k B T and the Fermi function is given by 



f(x) = 1/(1 + exp(x)). The contribution a q of the Dirac 
quasiparticles to the conductivity of a graphene mono- 
layer is then 
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where /j, is the chemical potential, and the factor of 4 
accounts for the spin and valley degrees of freedom. We 
find 
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where n is the scattering rate of the quasiparticles, which 
can be accounted for [16( by broadening the delta func- 
tions according to tt5 (x) = r]/(x 2 + rj 2 ). While the T- 
dependence of 77 is a priori unknown, we find that the 
scenario of constant /3rj, which may be ascribed to scat- 
tering off impurities or thermally generated ripples |16l | , 
is strongly favored by the available data in the range 
35 K < T < 150 K. The integral over £ in Eq. © can be 
performed analytically, and in the DC limit it yields 
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where z = e + in. The dependence of a q on /3A and fir] 
is illustrated in Fig. [1] 

We now analyze the data of Ref. Q on the suspended 
graphene devices SI, S2 and S3 in terms of the expression 
a = a + a b , where a b denotes the background con- 



ductivity. The simplest form compatible with the low- 
T data is the linear one a b = ct (1 — T /T). A more 
specific choice is the Variable-Range Hopping (VRH) 
law a b = a n expj — (T n /T) a ] encountered in conventional 
semiconductors [12| . While the former can be viewed as 
a linearized VRH expression, it cannot remain valid at 
arbitrarily low T as it becomes unphysical in that limit. 
As the T-dependence of the background is experimen- 
tally found to be weak, the choice of linear versus VRH 
description has very little impact on the analysis. 

In order to determine a b in an unbiased fashion, we 
first fix cr and T using the data in the low-T region 
where thermal activation is negligible. The next step is 
to subtract a b at all T, and determine j3r} and A by fit- 
ting a q to the resulting dataset. Finally, the results were 
confirmed by a simultaneous fit of all four parameters 
to the full dataset. In all cases negligible variation was 



observed, indicating that a Q and T are effectively uncor- 
rected with /3r] and A. The optimal parameter values 
are summarized in Table[Q and the results are plotted 
against the data of Ref. [6] in Fig. O 

Our findings indicate that the suspended graphene de- 
vices of Ref. [6| exhibit, upon subtraction of a b , a ther- 
mally activated component which can be well described 
in terms of Eq. © from T - 150 K down to T - 35 K, 
where the signal is lost due to limited measurement ac- 
curacy. As shown in Fig. [21 this corresponds to exponen- 
tial behavior over more than two orders of magnitude, 
with bandgaps in the range ^25 — 40 meV. The value 
j3r) ~ 0.1 obtained for samples SI and S2 is consistent 
with the high carrier mobilities and long mean free paths 
reported in Ref. [6]. Specifically, for T = 35 — 150 K we 
find 77 = 3.5—15 K, with corresponding mean free paths of 
hv F /r] ~ 2.0 — 0.5 /jm. For such low values of /3r], the rate 
of exponential decay is characterized by A whereas 0r) 
mainly determines the amplitude. It is noteworthy that 
the values of fir\ and A are roughly sample-independent. 
We have checked that fits with zero gap (A = 0), con- 
stant j] or zero background are incompatible with the 
data. Above T ~ 150 K, the data deviate from a de- 
scription with constant /3r), which may be ascribed to 
increasing phonon scattering at high T |6}. 

Once the description of the zero-density data is fixed, 
it is possible to predict a at finite carrier density n, with 
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where g ± (x,/3A,0ii) = f(^/x 2 + (/3A/2) 2 T Pn). It 
should be emphasized that o~(n) is strongly dependent 




FIG. 1: (Color online) DC conductivity a q as a function of /3A 
and f3n, shown semilogarithmically. For /3A S> 1, a q decreases 
exponentially when fin <g 1, while for fin S> 1 any signature 
of a gap is washed out by the large scattering rate. For fin <^. 
1, a„ decreases with fin, while for fin S> 1 this behavior is 
reversed. The shaded area denotes the approximate range in 
fin for suspended graphene samples. 



TABLE I: Optimal parameter values for the suspended 
graphene devices of Ref. Jy], corresponding to the analysis in 
Fig. [2] Similar results for the gap A and the scattering rate 
/3r] were obtained by subtracting the background conductiv- 
ity using a linear form or a VRH description with a ~ 1/3. 
Entries labeled by an asterisk (*) indicate that /3rj cannot be 
constrained due to lack of information on the normalization 
factor cr(5K), which is not known for device S3. The resulting 
uncertainty in A is negligible. 



sample 


Pv 


A[meV] 


crolkQ- 1 ] 


To [K] 


xV^dot 


Sl(a = l/3) 


0.103(3) 


36.8(1) 


0.530(6) 


0.027(7) 


1.6 


S2(a=l/3) 


0.105(1) 


26.2(1) 


0.485(5) 


0.85(6) 


1.1 


S3(a = l/3) 


0.070(1)* 


35.8(1) 


0.471(2) 


1.87(8) 


4.1 


S3(a = l/4) 


0.070(1)* 


36.4(1) 


0.551(2) 


2.81(8) 


4.0 


S3(linear) 


0.070(1)* 


32.0(1) 


0.339(1) 


2.21(1) 


7.1 



on 07], which is obtained via analysis of the quasiparti- 
cle contribution to cr(T) at n = 0. In Fig. |31 we show 
that the data of Ref. [6| at finite n is well described using 
the parameters of Table HI as obtained from the gapped 
quasiparticle analysis at n = 0. Notice also that the 
experimentally observed a ~ \/\n\ dependence is repro- 
duced (see Fig. [21 inset), as well as the transitional carrier 
density n* ~ 10 10 cm~ 2 which was identified in Ref. |6| as 
separating metallic (\n\ > n*) and insulating (|n| < n*) 
regimes in a(n,T). In the metallic regime, we recover 
the experimentally observed resistivity p ~ T. 

While a gap in the quasiparticle spectrum is an attrac- 
tive interpretation of the thermally excited conductivity, 
it should be noted that a number of other mechanisms 
can yield similar results. These include localized disor- 
der [l8|, scattering by screened charged impurities 19 1, 
and transport gaps due to quantum confinement 



11] 



Nevertheless, the gapped scenario is appealing as it pro- 
vides a straightforward and consistent description of the 
T- and n-dependence of cr, involving only quasiparticle 
excitations and an inhomogeneous regime close to the 
Dirac point. A plausible mechanism for gap generation 
in suspended graphene is given by the excitonic scenario, 
where a gap is dynamically generated by strong electron- 
electron interactions [20J,|21[. This idea has recently been 
revived based on Lattice Monte Carlo simulations of the 
low-energy effective field theory of graphene [22|, |23| . It 
is noteworthy that graphene flakes decoupled from un- 
derlying graphite layers have been found to exhibit a 
~ 10 meV bandgap at zero magnetic field, with a cor- 
responding splitting of the lowest Landau level of similar 
magnitude [8|. 

What is the origin of the observed background con- 
ductivity? Empirically, the physics at the Dirac point is 
obscured by charge density inhomogeneities [a, |7j , also re- 
ferred to as "puddles" , where n saturates to a finite value. 
The scale at which this happens for presently available 
samples is n ~ 10 11 cm -2 in non-suspended graphene, 
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FIG. 2: (Color online) Dirac Quasiparticle (QP) and back- 
ground (BG) components (see Table|I]) of the DC conductivity 
for the suspended graphene devices S1-S3, reproduced from 
Ref. [g]. All devices show a "knee" separating the thermally 
activated and background regions. After background subtrac- 
tion, device S3 exhibits exponential behavior over more than 
two orders of magnitude. The slight curvature at high T is 
due to the finite scattering rate fir\. 



and n ~ 10 9 cm 2 in suspended graphene [g, @|- One 
can define energy and temperature scales E and T via 



E 



k B T: 
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below which the description in terms of Dirac cones 
breaks down. Typically, E ~ 40 meV in graphene on 
a substrate, while in suspended samples E < 10 meV as 
h is reduced by an order of magnitude, thereby making 
it possible to access the physics at the Dirac point. Thus 
a bandgap of A ~ 30 meV is likely to be obscured by 
charge inhomogeneities in the non-suspended samples, 
whereas in suspended ones such a gap should be (par- 
tially) accessible. In samples with yet lower inhomogene- 
ity, one should observe a decrease in the low-T back- 
ground conductivity and an enhancement of the strongly 
T-dependent quasiparticle contribution. A description of 
<7 bg which is compatible with data is given by the VRH 
model 12], which describes the residual conductivity in 
semiconducting materials when thermal excitation is neg- 
ligible. If the background conductivity is indeed of the 
VRH form with a — 1/3 (as expected in 2D systems), 
it may indicate hopping between localized states. How- 




nen for instructive discussions and comments. Part of 
this work used the CUBPACK numerical quadrature rou- 
tine. 
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FIG. 3: (Color online) Left panel: 

dence of the DC conductivity a = a q + a bg for the sus- 
pended graphene device SI at T = 40 K for our model (solid 
red line) and experiment [g| before (triangles) and after (cir- 
cles) current- annealing. The dashed blue line corresponds to 
f}rj ~ 1.4, although the properties before annealing are likely 
dominated by charged impurities |6j. Inset: a as a function 
of \f n\. Right panel: Resistivity p(n,T) for device SI. As in 
Ref. @], dp/dT is metallic above n* ~ 10 10 cm -2 . 



ever, the data are inconclusive as the variation of a is 
very mild at low T. 

In summary, we have explored the signatures of a 
bandgap in the DC conductivity of graphene and showed 
that the empirical conductivity profiles a(n,T) are con- 
sistent with an interpretation in terms of Dirac quasi- 
particles with a non-zero bandgap A. However, the 
associated thermally activated behavior is partially ob- 
scured below an empirically observed characteristic scale 
h, where the carrier density as a function of the gate volt- 
age saturates. As gapped graphene is of great interest 
for nanoelectronic applications, a fully non-perturbative 
calculation (such as Lattice Monte Carlo) of the trans- 
port properties, including the effects of strong electron- 
electron interactions, is clearly called for. Such cal- 
culations also appear timely, as recent experimental 
work [l|, [2|, l24| at finite magnetic field has demonstrated 
a rich spectrum of phenomena closely related to Dirac 
physics at strong Coulomb coupling. 
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